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ABSTRACT 

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the 
EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydrody- 
namical cosmological simulations enables us to span four orders of magnitude in halo mass 
(11 ^ logio(Vf 2 oo/[ 1 M 0 ]) < 15), a wide radial range (—2.3 < logio(r/[/i _1 Mpc]) < 
1.3) and redshifts 0 z ^ 1. The shape parameters of the dark matter, stellar and hot gas dis¬ 
tributions follow qualitatively similar trends: they become more aspherical (and triaxial) with 
increasing halo mass, radius and redshift. We measure the misalignment of the baryonic com¬ 
ponents (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, 
redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align 
well with the local distribution of the total (mostly dark) matter. However, the stellar distribu¬ 
tions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to 
their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes 
(13 ^ logio(Af 2 oo/[/j _1 M 0 ]) ^ 15). Half of the disc galaxies in the EAGLE simulations 
have a misalignment angle with respect to their host haloes larger than 40 degrees. We present 
fitting functions and tabulated values for the probability distribution of galaxy-halo misalign¬ 
ment to enable a straightforward inclusion of our results into models of galaxy formations 
based on purely collisionless N-body simulations. 

Key words: cosmology: large-scale structure of the Universe, cosmology: theory, galaxies: 
haloes, galaxies: formation 


1 INTRODUCTION 

The topology of the matter distribution in the Universe is well de¬ 
scribed as a web-like structure comprising voids, sheets, filaments 
and haloes. This so-called cosmic web arises naturally from the 
gravitational growth of small initial perturbations in the density 
field of an expanding cold dark matter dominated (ACDM) Uni¬ 
verse. The evolution of the properties of the large-scale cosmic 
web is governed by the dominant components, i.e. dark energy and 
dark matter, while baryons are expected to trace the distribution 
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of the latter. Specifically, galaxies reside in dark matter haloes and 
trace them in terms of their positions and, to first order, in terms of 
their shapes and mutual alignment, albeit in a biased fashion due 
to the dissipative processes they experience during galaxy forma¬ 
tion. Theoretical stu dies of this galaxy bias have be en ongoing for 
several decades (e.g. lKaise ill 19841 : Da vis et al.ll 19851) . 

It has become apparent that when galaxies are used to infer 
the properties of the underlying dark matter distribution, it is con¬ 
venient to bisect this investigation into two steps: the relation be¬ 
tween galaxies and haloes and the relation between haloes and the 
underlying density field. The latter can be studied directly via cos¬ 
mological N-body simulations, whereas the former is a far more 
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complicated relation that is potentially affected by virtually all the 
physical processes associated with galaxy formation. For instance, 
while the triaxial shape of dark matter haloes is understood in terms 
of the collisionless nature of dark matter coupled with ellipsoidal 
collapse, galaxies manifest themselves in a plethora of morpholo¬ 
gies ranging from thin to bulge-dominated discs and to ellipsoidals 
and this is undoubtedly related to the redistribution of angular mo¬ 
mentum occurring during galaxy formation and evolution which, 
in turn, depends on the physical processes in operation. Thus, the 
characterization of the way galaxy shapes relate to their host haloes 
holds the potential to unveil the relevant physical mechanisms be¬ 
hind such a rich manifestation of galaxy types. 

Numerical simulations have been used to study the mutual 
alignment of g a laxi es with the© own host_ h aloes. For in stance , 
van den Bosch et all d2002j), IChen. Jing & Yoshikawl 120031). 


Sharma & Stein mew d2005li . iBett et al* d2010h and Sales et alj 
d2012h have shown that the angular momentum distributions of gas 
and dark matter components are partially aligned, with a typical 
misalignment angle of ~ 30°, although this might predominantly 
apply to disc galaxies. On the other hand, central ellipticals are 
expected to b e aligned with their host haloes i f they are f orme d 
by jnergers i Dubinski| 19981: Naab. Khochfar & Burkertl 1 20061 ; 


iBovlan-Kolchin, Ma & Ouataerl l^OOfil) . because the orientations 
of the central ellipticals and of the host dark matter haloes are 
determined by respectively the orbital angular momenta of their 
(correlated) progenitor galaxies and haloes. Observationally, there 
exist different indications of the presence of a misalignment 
between galaxies and their host haloes. Flowever, different studies 
have reached somewhat conflicting conclusion s abo ut the t ypical 
values of this m i salignment angle (see e.g. |Hgymans_eta]J 20041 : 


iKang et all2007l:lwang et all2QQ8l:lokumura. Jing & Lil2009l). 

Beyond its theoretical relevance, the misalignment of a galaxy 
with its own host halo can be a source of systematics for those stud¬ 
ies that aim to: infer the shape of dark matter haloes or constrain 
cosmological parameters via the measurement of the galaxy shape 
correlation function. Several current and forthcoming weak lensing 
surveys (e.g. KiDS, DES, LSST, and EuclicQ) will achieve the sta¬ 
tistical power to probe, observationally, halo shapes and to obtain 
exquisite measurements of the apparent alignment of galaxy shapes 
-cosmic shear- due to the gravitational lensing effect caused by the 
underlying (dark) matter distribution. It is therefore of great im¬ 
portance to guide the interpretation of the measured signal with 
numerical simulations. For instance, the link between the shape 
of the visible, baryonic matter and the structure of the underly¬ 
ing dark matter distribution, as well as their mutual orientation can 
be examined. To this end it is necessary to complement the ex¬ 
pectations derived from cosmological N-body simulations with the 
properties of galaxies as inferred from small-scale, high-resolution 


Joachimi et alj 20131: van den Bosch et al. 

2002 

: 1 Croft et alj 2009 

Hahn, Tevssier & Carollo 2010; Bett et al 

2010 

; Betl20li). 


In this paper, we extend previous work by exploiting the 
wealth of information encoded in hydro-cosmological simulations 
in which the main physical processes responsible for galaxy for¬ 
mation and evolution are simultaneously at play, thus leading to 


a more realistic realization of the galaxy-dark matter connection. 
We use the Overwhelmingly Large S i mulations (cosmo-OW LS, 

ISchave et alJl201(jl : lLe Bran et alJl2ol4l : lMcCarthv et alj|2014l) and 

the Evolu tion and Assembly of GaLaxies and t heir Environments 
(EAGLE. ISchave et alj|201a : ICrain et alj|2015t) project. This ap¬ 
proach has the advantage that the processes that lead to galaxy for¬ 
mation are self-consistently incorporated in the simulations and are 
therefore accounted for in the resulting galaxy and halo shapes, as 
well as in their correlation. During the late phase of this projec t, 
a study adopting a similar methodology was lTenneti et alj d2014l ). 
hereafter Tenl4, which has many aspects in common with our anal¬ 
ysis. Throughout the paper, we will therefore compare mutual find¬ 
ings. 

Our study is, however, unique as a consequence of sev¬ 
eral key features of our simulations and analysis. As detailed in 
§ 12.11 the use of cosmo-OWLS and EAGLE provides us suffi¬ 
cient cosmological volume and resolution, both of which are cru¬ 
cial for the reliability and the applicability of our results. Specif¬ 
ically, we span four orders of magnitude in halo mass (11 s( 
log(M 2 oo/[/i _1 M@]) L 15) and over six orders of magnitude 
in subhalo mass M su b, enabling us to investigate spatial variations 
of the shape of galaxies and haloes from galactic to cosmologi¬ 
cal scales. Furthermore, the combination of EAGLE and cosmo- 
OWLS forms a set of simulations that reproduces the observed 
abundance of galaxies as a function of stellar mass (the galaxy stel¬ 
lar mass function) at both low (log(M 2 oo/[ft 1 Mg]) ^ 13) and 
high (13 < log([M 2 oo/[fi 1 M@]) ^ 15) halo masses. Moreover, 
it has been shown that the cosmo-OWLS simulations reproduce 
various (X-ray and optical) observed properties of galaxy gro ups 
dCrain et al.1 l20ld : iMcCarthv et alj 1201 d : iLe Bran et al.112014l) as 
well as the observed galaxy mass function for haloes more massive 
than log(M 2 oo/[ (i _1 Mg]) = 13. Finally, the gal axy size distribu¬ 
tion in EAGLE reproduces the observed one (ISchave et all2015h . 

This paper is organized as follows. We summarize the proper¬ 
ties of the simulations in §[2] where we also introduce the technical 
definitions used throughout the paper. We highlight some caveats 
to the shape and angle estimates related to the feedback implemen¬ 
tation in §[3] In §[4]we present the results concerning the sphericity 
and triaxiality of dark matter haloes, as well as those of the stellar 
and the hot X-ray emitting gas distribution. The (mis)alignment of 
the baryonic components with their host haloes is addressed in §[5] 
We summarize and comment on our results in § [6] 

Throughout the paper, we assume a flat ACDM cosmology 
with massless neutrinos. Such a cosmological model is character¬ 
ized by fiv^| parameters: {D m , fib, erg, n s , h}. The simulations 
used in this paper were run with two slightly different sets of values 
for these parameters. Specifically, we will refer to PLANCK as the 
set of cosmological values suggested by the Planck mission (D m , 
D b,08, n s , h } = 10 .307, 0.04825 0.8288, 0.9611, 0.6777} (Table 
9: IPlanck Collaboration et alj2014l) . whereas WMAP7 refers to the 
cosmological parameters (fl m , fib, &8, n s , hj = (0.272, 0.0455, 
0.728, 0.81, 0.967, 0.704} suggested by the 7th-year data release 
dKomatsu et alj20lH) of the WMAP mission. 


1 KIDS: Kilo-Degree Survey, 
http://www.astro-wise.org/projects/KIDS/ 

DES: Dark Energy Survey, https://www.darkenergysurvey.org; 

LSST: Large Synoptic Survey Telescope, http://www.lsst.org 

Euclid: http://www.euclid-ec.org 2 Flatness implies that Q,\ = 1 — Q m . 
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Table 1. Simulations used throughout the paper and their relevant properties. Description of the columns: (1) descriptive simulation name; (2) co¬ 
moving size of the simulation volume; (3) total number of particles; (4) cosmological parameters used in the simulation; (5) initial mass of baryonic 
particles; (6) mass of dark matter particles; (7) maximum softening length; (8) colour used for the simulation; (9) simulation name tag. 


Simulation 

(1) 

L 

(2) 

-^particle 

(3) 

Cosmology 

(4) 

m h 

[ft- 1 M 0 ] 
(5) 

m dm 

[h- 1 Mq] 
(6) 

Cprop 

(ft -1 kpc) 

(7) 

Colour 

(8) 

tag 

(9) 

EAGLE Recal 

25 (Mpc) 

2 x 752 3 

PLANCK 

1.5 x 10 s 

8.2 x 10 5 

0.5 

purple 

EA L025 

EAGLE Ref 

100 (Mpc) 

2 x 1504 3 

PLANCK 

1.2 x 10 6 

6.6 x 10® 

0.2 

orange a 

EAL100 

cosmo-OWLS AGN 8.0 

200 (ft” 1 Mpc) 

2 x 1024 3 

WMAP7 

8.7 x 10 7 

4.1 x 10 s 

2.0 

blue 

CO L200 

cosmo-OWLS AGN 8.0 

400 (ft” 1 Mpc) 

2 x 1024 3 

WMAP7 

7.5 x 10 8 

3.7 x 10 9 

4.0 

green 

CO L400 


a Cyan is used for Figs. II ll and ll2l where the EA L100 simulation is used in order to improve the statistics for the least massive bin. 


2 SIMULATIONS AND TECHNICAL DEFINITIONS 


2.1 Simulations 

Throughout the paper, we employ the outputs of four cosmologi¬ 
cal volu mes simulated within the context of tw o distinct projects: 


EAGLE l Schaye_etaL 2015^ Crain et akll2015l) and cosmo-OWLS 
jLe Bran et aljl2014l : lMcCarthv et al.ll2014h . We use the former to 
investigate (well-resolved) smaller halo masses in relatively small 
volumes; whereas the latter is used to study more massive haloes in 
larger volumes. Table Q] lists all relevant specifics of these simula¬ 
tions. 

Both EAGLE and cosmo-OWLS were run using a mod¬ 
ified version of the iV-Body Tree-PM smoothed particle hy¬ 
drodynamics_(SPH) code GADGET 3, which was last de¬ 

scribed in ISpringell J2005h . The main modifications are the 
formulation of the hydrodynamics, the time stepping and, 
most importantly, the subgrid physics. All the simulations 
used in this work includ e ele ment-b y-ele ment radi ative cool¬ 
ing (for 11 eleme nts: Iwiersma. Schave & Smith] l2009h . star 
formation i Schave & Dalla Vecchial 120081) . stellar mass loss 
l Wiersma_et_aL 2009|), energy feedb ack from star formation 
l Dalla Vecchia & Schavel 120081. 120121) . gas accretion_ onto and 
mergers of supermassive b lack holes (BHs: iBooth & L SchaY j2009|: 
Ros^-GuevarartaL 20131), and AGN feedback ( Booth & Schavel 
20091: ISchave et al.ll2015l). 

The subgrid physics u sed in EA GLE builds upon th at of 
OWLS dSchave et all l20ld) GIM I C dCrain et all l2009al) and 
cosmo-OWLS dLe Bran etall 12014 iMcCarthv etaU 12014 . Fur¬ 
thermore, the EAGLE project brings a number of changes with 
respect to cosmo-OWLS regarding the implementations of energy 
feedback from star formation (which is now thermal rather than ki¬ 
netic), the accretion of gas onto BHs (which now accounts for an¬ 
gular momentum), and the star formation law (which now depends 
on metallicity). More information regarding technical implementa¬ 
tion of hydro-dy namical asp ects as well as subgrid physics can be 
found in lschave et al.l J2015I) . 

Arguably, the most important feature of the EAGLE simu¬ 
lation is the calibration of the subgrid physics parameters to re¬ 
produce the observed galaxy mass function and galaxy sizes at 
redshift zero. One of the key feature of the cosmo-OWLS simu¬ 
lations is that they reproduce optical and X-ray scaling relations 
of groups and clusters of galaxies. In this work we exploit both 
these unique features by splitting our range of halo masses into four 
mass bins and by using a different simulation for each one of them. 
Specifically, for halo masses below the ‘knee' of the galaxy stel¬ 
lar mass function we use EAGLE in order to ensure galaxies form 


with the ‘correct’ efficiency and size, whereas for haloes above 
the ‘knee" we use cosmo-OWLS. In practice, we create a compos¬ 
ite sample of haloes spanning four orders of magnitude in mass 
(11 < log(M 20 o/[/i _1 Mq]) < 15). 


2.2 Halo and subhalo definition 


Groups of particles are identified in our simulations by applying 
the Friends-of-Friends_ algorithm with linking length 0.2 to the dark 
matter particles J Davis et aD 19851) . The mass M 2 oo and the radius 
r 200 of the groups are assigned using a spherical over-density al¬ 
gorithm centred on the mi nimum of the gravitational potential, a s 
implemented in SUBFIND JSnrinsiel et al.1200ll ; lDolag et alj|2009h . 
From each group, dynamically un-bounded particles are discarded. 
Thus, subhaloes are identified as a collection of bound particles that 
reside in a local minimum of the gravitational potential computed 
using all particle types. The most massive subhalo is the central 
subhalo of a given FoF group and all other subhaloes are satellites. 
Particles that are bound to a subhalo belong exclusively to that sub¬ 
halo. Correspondingly, central subhaloes do not contain particles 
that reside in other local minima of the potential, even if those par¬ 
ticles are within the subhalo boundary. We define the centre of a 
subhalo as the position of the particle with the minimal gravita¬ 
tional potential. The subhalo radius can be calculated for each com¬ 
ponent separately. A commonly used estimate is the radius within 
which half of the mass in dark matter is included, rjJSf- The mass 
of a subhalo is the sum of the masses of all the particles that consti¬ 
tute it. For the rest of the paper we will use the term ‘halo’ to refer 
both to central and satellite subhaloes, unless otherwise specified. 

The masses of subhaloes for both centrals and satellites (ac¬ 
cording to SUBFIND classifications), are indicated with M su b. 
However, whenever a distinction is required, we shall use M 200 
and r 2 oo to characterize the properties of central haloes. 


2.3 Shape parameter definitions 

A fundamental quantity that describes how matter is spatially dis- 
tributed is the three-dimensional ma ss distribution tensor ( e.g. 
Da vis et anil985l : ICoie & Lacevll 19961) . 

iVp 

Afjj — ^ ^ rripXpiXpj , (1) 

p=i 

where Np is the number of all particles that belong to the structure 
of interest, x p i denotes the element i (with i,j = 1,2,3 for a 3D 
particle distribution) of the position vector of particle p, and m p is 
the mass of the p th particle. This mass distribution tensor is often 
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Table 2. Values at z = 0 of various quantities of interest in each mass bin. Description of the columns: (1) simulation tag; (2) mass range 
logio (M 20 Q /[h~i M©]) of the haloes selected from the simulation; (3) median value of the halo mass logiofM^QQ); (4) median value of the 
stellar mass (logio(M s tar/[ h~ 1 Mq])) considering all the star particles that belong to the halo; (5) standard deviation of the stellar mass distribution 
^logioMstar» (6) median value of halo radius r^jjo; 0) median radius within which half of the mass in dark matter is enclosed; (8) median radius 
within which half of the mass in stars is enclosed; (9) number of haloes; (10) number of satellite haloes. 


Simulation tag 

mass bin 

* 

mho 

* 

-Mstar 

* 

^logio-Mstar 

* 

^crit 

'200 

** 

^dm 
' half 

** 

~,star 
' half 
** 

-^Yhalo 

N s &t 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

EA L025 

[11 - 12] 

11.31 

9.50 

0.45 

96.0 

39.8 

2.7 

156 

24 

EA LI 00 

[12 - 13] 

12.27 

10.58 

0.26 

199.3 

93.4 

4.9 

1008 

104 

CO L200 

[13 - 14] 

13.16 

11.21 

0.25 

396.4 

241.8 

53.4 

2190 

137 

CO L400 

[14 - 15] 

14.09 

12.06 

0.19 

805.9 

505.1 

106.7 

1152 

26 


* logio(M/[ ft 1 M 0 ]) 
** 77/[ft _1 kpc] 


referred as the in ertia tensor, since the two tensors share the same 
eigenvectors (see lZemp et alj|201 lL for a discussion) and, for most 
astrophysical purposes, those eigenvectors encode the information 
of interest. Throughout this paper we will refer to the mass distri¬ 
bution tensor as the inertia tensor to conform to the jargon used in 
the literature. 

The eigenvalues of the inertia tensor will be denoted as Ai 
(with i = 1,2,3 for a 3D particle distribution as in our case). Given 
a particle distribution inertia tensor, the modulus of the major, in¬ 
termediate, and minor axes of the corresponding ellipsoid can be 
written in terms of these eigenvalues as a = y/Ay, b — y/\ 2 , and 
c = vA 3 > respectively. We interpret this ellipsoid as an approx¬ 
imation to the shape of the halo. Specifically, the sphericity and 
triaxiality parameters, S and T, are defined as 

S= and T = ° 2 ~ b l . (2) 

a or — c A 

A purely spherical halo will have S = 1 with T being undefined. 
Low values of T (i.e. T — > 0) correspond to oblate haloes while 
high values (i.e. T — > 1) correspond to prolate haloes. 

We note that the computation of shape parameters in a spher¬ 
ical region biases the shape towards higher sphericity. When com¬ 
puting shapes of dark matter haloes in spherical regions it is pos¬ 
sible to correct for this effect applying the simple empirical re¬ 
scaling: Strue = S^ as suggested in lBailin & Steinmetzl d2005h . 
This correction is not implemented in the results presented here 
since a similar correction is not available for the other quantities 
that we present. 

2.4 Axes and misalignment angle definition 

The eigenvectors of the inertia tensor, in Eq[T| are denoted as ej, 
with i = 1, 2,3 in the case of a 3D distribution of particles and 
x = halo, star, gas to indicate total matted stars, or gas, respec¬ 
tively. We relate the ordered eigenvectors ej~, e^ 2 , and e., 3 of the 
inertia tensor to the direction of the major, intermediate and minor 
axis of the corresponding ellipsoid. We further indicate the radial 
dependence of the major axis as e\(r), which, unless stated other¬ 
wise, has been computed using the volume enclosed by the entire 
structure as defined by SUBFIND (see 112.21 . We shall quantify the 

3 We do not deal with the specific case of only dark matter because on the 
scales of interest it almost exactly coincides with the total matter in a halo. 



l° 9 io ^200 [M 0 ] 

Figure 1. The stellar mass to halo mass ratio of central galaxies as a func¬ 
tion of halo mass, normalised by the cosmic baryon fraction, for the four 
feedback variations used in J3] The curves are dotted where there are fewer 
than 100 star particles per galaxy and individual galaxies are showed for 
bins that contain fewer than 10 galaxies. The 1-sigma scatter about the me¬ 
dian of Reference is shown as a shaded region. Dark and light grey l ines rep¬ 
resent the abundance matching relations of lBehroozi, Wechsler & Conrovl 
l2013l f and Moster, Naa b & White |20 1 3| ). 

alignment of different matter components via the scalar product of 
two major axes, i.e. the misalignment angle 9 (0 in case of pro¬ 
jected quantities). Specifically, we will use cos 9 as the principal 
quantity of interest and only comment on the actual value of 9 when 
relevant. We stress here that the major axis is a spin-2 quantity, i.e. 
it is invariant under rotation of 180 degrees. This means that 9 only 
varies between 0 and 90 degrees and, correspondingly, cos 9 can 
only assume values between zero and unity. 


3 THE EFFECT OF GALAXY FORMATION 
EFFICIENCY 

A major asset of our composite sample of simulated haloes is that 
it reproduces the observed stellar-to-halo mass ratio as a function 
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Figure 2. Ratios of the average sphericity (left panel) and average misalignment angles (right panel) of the stellar distribution with respect to the average 
values of the REFERENCE EAGLE simulation (see m . Different colours indicate different mass bins while different line styles refer to different simulations 
which differ only by the implementation of feedback. The disagreement (up to 20% for sphericity and up to 40% for misalignment angle) stems from the 
different efficiency of galaxy formation (see discussion in §[3j. The misalignment angle is more sensitive to galaxy formation efficiencies than the sphericity. 
The differences always increase towards the centre of the halo. 


of halo mass. Specifically, EAGLE has been calibrated to repro¬ 
duce the stellar mass function at redshift zero and cosmo-OWLS 
has proven successful in reproducing many observable propertie s 
of groups and clusters ^McCarthy et al.ll2Qlol : lLe Brun et all2014h . 
Moreover, in the halo mass range where cosmo-OWLS haloes are 
us ed, th eir .galaxy formatio n effic iency is consistent with the results 
of Moster. Naab & Whitel d2013l) : iBehroozi. Wechsler & Conrovl 
( 2013 ) from abundance matching techniques. This feature is par¬ 
ticularly important in the context of our investigation, as one might 
expect that if a simulation produces either too many or too few 
stars, then their distribution and consequently, the galaxy shape pa- 
ra meters w ould also be affected. Note however that, as shown by 
I Crain et aD ( l2015h . this criteria is insufficient to guarantee that the 
spatial distribution of barionic matter is realistic. 

Before showing the main results of our analysis we investi¬ 
gate how different feedback implementations results in different 
predictions for the shape and orientations of galaxies with respect 
to their host haloes. To quantify this effect, we make use of a set 
of feedback variations on the Reference model of the EAGLE sim¬ 
ulations that, unlike the Reference model itself, do not reproduce 
the observed galaxy stellar mass function (i.e. in those simulations 
haloes do not form stars with the observed efficiency). A detaile d 
description of these simulations can be found in lCrain et alj d2015h . 
Here, we only briefly summarize their properties. All simulations 
adopt the PLANCK cosmology. The simulation boxes have comov¬ 
ing volumes of 25 3 Mpc 3 , with 2x376 3 particles. We consider four 
variations: 


L025_ref A simulation run in a smaller volume with respect to the 
main run using the Reference EAGLE implementations, namely: 
thermal energy feedback from star formation, BH gas accretion that 
takes into account the gas angular momentum and a star formation 
law which depends on gas pressure and metallicity. In the thermal 
feedback implementation the amount of energy injected per feed¬ 
back event is fixed but there is freedom in the amount of energy 
that can be injected per unit of stellar mass. This freedom is in¬ 


corporated in the parameter fth that is the expectation value of the 
amount of energy injected per unit stellar mass formed, in units 
of the energy available from core collapse supernova for our IMF. 
The average number of neighbouring particles heated by a feedback 
event is < Wheat >~ 1.3/th ( io j .Ty ) 1 whereas the temperature 
jump for the single particle is fixed to AT = 10' ° K. If the value 
of fth is constant, then both the energy injected per single event 
of feedback and the energy per unit of stellar mass are fixed. By 
varying the parameter fth, it is possible to control the efficiency of 
the feedback and so to account for the unresolved radiative losses 
that depend on the p hysical state of the IS M, or to compensat e for 
numerical losses (see lSchave et al . 20H and Crain eTal l2015l for a 
discussion). The value of fth depends on the local physical condi¬ 
tions (density and metallicity) of the gas according to: 


fth — /th,min T“ 


fth ,max — fth ,n 


1 I ( Z Y Z / "H,birth V 

+ V°1 Z 0/ V n H,0 ) 


(3) 


where jih, birth is the density of the parent gas particle at the time it 
was converted into a stellar particle and Z is the gas metallicity. The 
value of riH,o = 0.67 cm -3 was chosen to reproduce the observed 
present-day GSMF and galaxy sizes, whereas nz = n n = 2/lnlO. 
We use the asymptotic values /th.max = 3 and /th.min = 0.3, 
where the high asymptote /th,max is reached at low metallicity and 
high density. 

L025_wfb Weaker stellar feedback than for the Reference model. 
In this case the function in Eq.[3]is scaled by a factor of 0.5. 
L025_sfb Stronger stellar feedback than for the Reference model. 
In this case the function in Eq.[3]is scaled by a factor of 2. 
L025_nag Same as Reference but without AGN feedback. 


Fig-0shows the stellar mass to halo mass ratio of central galaxies 
as a function of halo mass, normalized by the cosmic baryon frac¬ 
tion, for the four aforementioned feedback variations. The galaxy 
stellar mass function and the galaxy sizes as obtained from these 
different feedback variations can be seen in Fig. 10, panels a and 
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c, of ICrain et al .1 d2015h . Those models produce stellar mass func¬ 
tions with differences of the order of half a dex above (L025_wfb) 
and below (L025_sfb) the Reference one. The case without AGN 
feedback differs from the Reference case only for the most mas¬ 
sive galaxies. Dark and li ght gre y line s represe nt the abundance 
matching relation s of iBehroozi. Wechsler & Conrovi d2013h and 
iMoster. Naab & Whitel d2013l) . respectively. The Reference simu¬ 
lation shows good agreement with the abundance matching mod¬ 
els. 

Fig. [2] shows the changes in the main quantities of interest in 
our analysis for the aforementioned feedback implementations. The 
left panel displays the ratio of the sphericity of the stellar compo¬ 
nent of haloes as a function of the distance from the halo centre 
for each simulation with respect to L025_ref. Different line styles 
refer to different simulations and we report the results for two halo 
mass bins. The difference are of the order of 10%. For the triaxial- 
ity parameter (not shown here) the differences range from 15% to 
50 %. The right panel displays the ratio of the cosine of 6(r) (the 
angle between the halo’s first eigenvector and the first eigenvector 
of the stars inside a given radius) of each simulation with respect 
to L025_ref. This quantity shows 10% differences at r^oo, while 
differences as large as 40% for the case without AGN (and 20% in 
the case of weak SN feedback) are present on scales representative 
of typical galaxy sizes. We report that the differences between the 
sphericity of haloes in the different sub-grid implementations (not 
shown) are smaller than 5% at all radii. This analysis underlines the 
importance of the calibration of feedback, especially for the shape 
and alignment of the innermost parts of haloes where most of the 
stars reside. 

A priori, there is no guarantee that reproducing the galaxy 
stellar mass function is a sufficient condition to predict realistic 
shape parameters. For instance, one may envision a scenario in 
which the si ze of g alaxi es, at the same mass, will also influence 
their shapes. ICrain et ail ( 120151 ) have reported four different sim¬ 
ulations in which the galaxy stellar mass function is equally well 
reproduced but the predictions for galaxy sizes are widely differ¬ 
ent. We computed the shape parameters and star- halo misalignment 
for the same simulations employed in lCrain et ali J 2 OI 5 I) . Although 
in rough agreement, the relative variance from model to model is 
~ 10% — 15% for both the sphericity and the misalignment angle 
(not shown). Clearly, beyond the effect of the ‘galaxy formation 
efficiency’, galaxy sizes also play a role in the accuracy of the re¬ 
trieved shape parameters. 

In this section and in the rest of this paper we will not focus on 
the origin of the different shapes and misalignment of the different 
populations of haloes. Investigating the physical origin of shapes 
and misalignments represents an interesting line of inquiry that has 
been addressed using zoom-in simulations and by foll owin g the 
evolu t ion of galaxies and haloes in time ( e.g. iRomano-Dfaz et al] 
120091 : IScannanieco et al]|2009t ICenll2014h . In this work we focus 
on exploiting the large dynamical range available to give statistical 
trends with halo mass and radius and postpone a detailed investiga¬ 
tion on their physical origin to future work. 


4 SHAPE OF THE DIFFERENT COMPONENTS OF 
HALOES 

Armed with the simulations described in § 12.11 and with the tech¬ 
nical definitions introduced in §2.2-2.4, we present here a system¬ 
atic study of the shape parameters. Specifically, we will present 
the shape parameters of the entire matter distribution in haloes (in 


m. and of different halo components (stars in 1 14.2l and hot gas 
in i j4.3l l as well as their mass, spatial, and redshift dependence. In 
Table [3] we summarize the values and the scatter of the shape pa¬ 
rameters and the misalignment angles for the stars within . 

It is well known that the reliability of shape estimates of par¬ 
ticle distributions depends on the number of particles used to trace 
those distributions (e.g. Tenl4). Motivated by the results presented 
in Appendix lA2l we measure shapes of structures with at least 300 
particles. The resolution criterion is applied separately to the differ¬ 
ent halo components. Therefore, for a reliable shape measurement 
of the stellar component we require galaxies containing at least 300 
stellar particles. Our tests performed using synthetic NFW haloes 
show that this choice ensures a precision of 3% and an accuracy 
better than 10% in the sphericity and triaxiality parameters, see Ap¬ 
pendix lA2l for more details. We note that Tenl4 performed a sim¬ 
ilar convergence test according to which using 300 particles leads 
to ~ —10% bias in the sphericity of a particle distribution. Our 
choice ensures relatively high precision while still allowing us to 
have a large number of haloes for which shape measurements can 
be performed. 


4.1 The shape of haloes 

Rg. a displays the sphericity (left panel) and triaxiality (right 
panel), S and T respectively, for halo masses in the range 9 -s) 
logio(M su b/[hA 1 Mpc]) ^ 15. Different colours indicate differ¬ 
ent simulations and different line styles represent different redshifts 
(see legend). Notably, despite their difference in resolution, the re¬ 
sults agree in the overlapping mass intervals probed via different 
simulations. The common qualitative result is very simple: haloes 
become less spherical and more triaxial (prolate) with increasing 
mass. Sphericity (triaxiality) decreases (increases) from z = 0 
(solid lines) to z = 1 (dotted lines). Haloes thus become more 
spherical/oblate as cosmic time progresses. This effect is not due to 
baryon ph ysics since it was also found in dark m atter only simula¬ 
tions (e.g. lBrvan et al.ll2013l : lTenneti et alj|2014h. For comparison , 
we also plot the halo sphericity reported by iTenneti et al] d2014l) 
using a dashed line for z = 0 and a long dashed line for z = 1. 
Despite the differences in box size, resolution and implementation 
of baryon physics, the overall agreement with our composite set of 
simulations is good at both redshifts. The shape of the haloes when 
all particles are considered is dominated by the dark matter com¬ 
ponent. In fact, the shape of the dark matter component is nearly 
identical to the that of the total mass distribution. 

Our composite sample suggests that, over a wider range in 
halo masses, the relation deviates from linear showing a steepening 
from low to high masses. 

4.2 Shape of the stellar component of haloes 

Fig. 0 displays the halo mass dependence of the shape parameters 
of the stellar distributions. Sphericity is on the left, triaxiality is on 
the right. As in Fig. [3] different colours indicate different simula¬ 
tions according to TableQ] We remind the reader that we use a min¬ 
imum of 300 particles to determine the shape of particle distribu¬ 
tions. This inevitably leads to a relatively small halo mass range for 
each simulation. However, the composite sample of our simulations 
covers the halo mass range 11 ^ logio(M su b/[/i^ 1 Mg]) ^ 15. 
Note that we have indicated with grey lines the values of the 
shape parameters obtained when considering haloes comprising 
fewer than 300 particles. Interestingly, in the overlapping halo 
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Table 3. Values of main quantities of interest for each halo mass bin. Values refer to z = 0 and are measured at the half-mass radius in star, for 
all subhaloes. Angle 6 refers to 3D quantities, whereas © refers to the 2D projected equivalent. Description of the columns:(l) simulation tag; (2) mass 
range of the haloes, logio M 200 ; (3) Median value of the subhalo mass, logio M su b-> considering the sum of all the masses of the particles belonging 
to the subhalo; (4) median value of the stellar mass considering all the star particles belonging to the halo; (5) median value of the sphericity computed 
at the stellar half-mass radius; (6) median value of the triaxiality computed at the stellar half mass radius; (7) median value of the projected ellipticity 
(averaged over the three axis projections x, y and z); (8) median angle between the first eigenvector of the stellar component enclosed in and the 
first eigenvector of the total matter distribution in the halo; (9) same as (8) but for the projected haloes averaged over the three projection axes; (10) 
median angle between the first eigenvector of the stellar distribution and the total matter distribution, both evaluated at (11) same as (10) but for 
the projected haloes averaged over the three projection axes. 
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Figure 3. Halo sphericity (left) and triaxiality (right) as a function of halo mass. Both central and satellite haloes are considered, hence the choice of M su b 
(the sum of the masses all the particles belonging to the subhalo) as identifier of the halo mass. Both shape parameters are computed using all the particles 
in the subhaloes (gas, stars, and dark matter). Using only dark matter would give virtually identical results. Different colours indicate different simulations, 
whereas solid (dotted) lines refer to z = 0 (z = 1). The error bars represent the one sigma bootstrap error on the median. Dashed black lines are the values 
obtained using the fitting functions from Ten 14. 


mass range, the sphericity parameters derived from simulations 
with different resolutions agree remarkably well. The general trend 
seems to suggest that sphericity is a decreasing function of halo 
mass for logio(M su b/[ hT 1 M©]) > 12 at z = 0 and for 
logio(M sub /[/i _1 M 0 ]) > 11 at 2 = 1. 

We compa re our results in Fig. [4] with the recent work 
of iTenneti et al.l | |2014| ) by showing their fitting function to the 
sphericity of the stellar component of haloes (black dashed line). 
The most prominent feature of their fitting function, namely the 
sharp upturn at masses logio (M sub /[/i _1 M@]) < 11, is most 
likely due to a selection bias. In their work they only compute 
shapes for subhaloes with more than 1000 stellar particles. This 
choice imposes a strict limit in stellar mass but not in subhalo mass. 


This approach only results in an unbiased selection if the mini¬ 
mum stellar mass of all haloes in a given mass bin is higher than 
> 1000m s tar where m s tar is the mass of a stellar particle. If we 
impose the same strict limit of 1000 star particles without also lim¬ 
iting the halo masses accordingly, we obtain a similar upturn in the 
stellar sphericity. Moreover, this upturn occurs at a different mass 
for different simulations since a fixed number of particles translates 
into different mass depending on the resolution used. 

The triaxiality parameter (right panel of Fig. 0 is an increas¬ 
ing function of halo mass at both z = 0 and 1. As discussed in 
Appendix IA2I the accuracy of the triaxiality estimate is more sen¬ 
sitive to the minimum number of particles used to compute it. This 
manifests itself in the fact that the grey lines in this plot do not con- 
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Figure 4. Stellar shape parameters (sphericity on the left, triaxiality on the right) as a function of halo mass. Different colours indicate different simulations, 
whereas different line styles refer to different redshifts. The error bars represent the one sigma bootstrap error on the median. Grey lines are show the results 
for mass bins containing haloes with less than 300 stellar particles. The dashed black line indicates the sphericity obtained from the fitting function of Tenl4. 
The upturn and the downturn in this fitting function are likely due to selection effects (see discussion in § 14.2.21 . 



Figure 5. Projected rms stellar ellipticity as a function of halo mass. Dif¬ 
ferent colours indicate different simulations (and therefore halo masses), 
whereas different line styles indicate the region within which the stellar dis¬ 
tribution is considered. Specifically, dashed lines indicate the case in which 
only star particles within the entire halo are considered, whereas solid lines 
indicate the case in which star particles within the stellar half-mass radius 
are considered. Both centrals and satellites are considered this analysis. 


tinue a monotonic trend beyond the well-resolved mass interval, 
thus reinforcing the importance of imposing a minimum number of 
particles used when attempting to recover the triaxiality of a distri¬ 
bution of particles. 


4.2.1 The projected stellar rms ellipticity 

Under the assumption that galaxies are randomly oriented, aver¬ 
aging the observed projected ellipticity of galaxies gives a mea¬ 
surement of the gravitational lensing effect that, in turn, gives 
constraints on the matter distribution along the line of sight. 
The S/N of those measurements depend on the second mo¬ 
ment of the distribution of galaxy intrinsic ellipticity, termed 
Crms • Many observational studies have measured the value of 
the e rm s f or populations of gala xies . Early results were reported 
in iHoekstra. Franx & Kuiikeij hood , and more stati stically ro¬ 
bust result s were obtained using SPSS data dRev es et al.l l2012l) . 
COSMOS jjoachimietaj [201j : ^^ndelbaum^alj|2014 ) and the 
CFHTLenS survey i Hevmans et alj 2013 : Miller et all 20131) . Un¬ 
fortunately, despite the tremendous progress in the statistical power 
of the galaxy surveys employed in these studies, obtaining an ac¬ 
curate estimate of e rm s remains challenging, especially because 
of the fact that the quantity that is accessible observationally 
al ways has a (often non-negligible) noise contribution (see e.g. 
IViola. Kitching & Joachimill2014l) . 

For our composite sample of haloes, e r ms is defined as: 


= iV 

N ^ 


1 - <(C 

T+gf 


(4) 


where q[ is the projected ellipticity of the i th halo q' = h'/a' where 
a' and b' are the values of the major and minor axis of the projected 
stellar distribution, and N is the total number of haloes considered. 

We use our composite sample of haloes to compute the stellar 
firms in bins of halo mass of width 0.5 dex, as a function of halo 
mass in Fig [5] We make use of all star particles that belong to the 
subhaloes (dashed lines) or only stellar particles within the stellar 
half-mass radius (solid lines). Both centrals and satellites are con¬ 
sidered for this analysis. When all star particles are considered the 
value of the e rm s increases with mass from 0.35 to 0.55. System¬ 
atically lower values are found if only stellar particles within the 
half-mass radius are considered. 
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The values of the e rms predicted by our composite sample, 
when all stars are considered, are in broad agreement with the ob¬ 
served noise-corrected values that are of the orde r of sa 0.5-0.6 de - 
pending on luminosity and galaxy type (e.g. ljoachimi et aPl2013l) . 
Unfortunately, a direct comparison of our results with those ob¬ 
tained from observational studies is far from trivial. In fact, it would 
be crucial to mimic all steps in the observational methodology. 
For instance, e rms measurements are usually only available for a 
given sub-population of galaxies, those galaxies are further binned 
in absolute magnitude, and the axis ratio is computed starting from 
(noisy) images for which flux isophotes need to be identified. In the 
context of this investigation, we find the current level of agreement 
satisfactory and ideal as a starting point for future explorations. 


4.2.2 Variation of the shape of the stellar component of haloes 
with the distance from the halo centre 

Fig. [6] shows the sphericity (left panel) and the triaxiality (right 
panel) of the stellar component of haloes as a function of the dis¬ 
tance from the centre of the halo. 

We divide our sample into mass bins that are drawn from dif¬ 
ferent simulations according to Table |T| We then compute the iner¬ 
tia tensor for increasingly larger spheres around the centre of each 
halo. For every sphere we show the median values of the shape pa¬ 
rameters of the mass inside the sphere. Radii are given in units of 
r^oo to allow for a comparison of haloes of different masses. Only 
particles that are bound to the halo are considered for this analysis. 
Curves are drawn only on scales where at least 300 particles can be 
used. 

The stellar component of haloes tends to be more spherical 
near the centre. The triaxiality value shows significant evolution 
for masses below M su b < 10 12 hr 1 Mg. These trends are qual¬ 
itatively the same as those found for the dark matter component 
(not shown) with the exception that the radial profile of the stellar 
distribution is steeper than that of the dark matter distribution. 

The right panel of Fig. [6] shows a large difference be¬ 
tween the triaxiality values of subhaloes in the mass bins 12 < 
logio(M su b/[/i^ 1 Mq]) < 13 (orange curves, EA L100) and 
13 < logio(M su b/[ hr 1 M 0 ]) < 14 (blue curves, CO L200). 
This feature might be caused by the different resolution, vol¬ 
ume, and/or baryon physics of the two sets of simulations (al¬ 
though the latter is relatively small). To test whether that is the 
case, we compute the triaxiality parameter of subhaloes with 
mass 13 < logio(M au b/[/i _1 Mg]) < 14 using the EAGLE 
L100 simulation (not shown). We find the corresponding re¬ 
sults to agree with the results obtained using the same mass 
bin from the cosmo-OWLS L200 simulation. Thus, we inter¬ 
pret the differences between the triaxiality of subhaloes in the 
mass bins 12 < logio(Af su b/[ /i _1 Mg]) < 13 and 13 < 
logio(M su b/[/i^ 1 Mg]) < 14 as having a physical origin rather 
than being due to the resolution, the volume, or the (small) differ¬ 
ences in the baryon physics of the two sets of simulations. 

4.3 Shape of the hot gas component of haloes 

In this section we repeat, for the hot gaseous component of haloes, 
the analysis performed for the total and stellar matter in § 14.11 and 
§ 14.21 respectively. We present the shape parameters for a subsam¬ 
ple of temperature-selected diffuse gas (T > 10 6 K). The selection 
is quite insensitive on the exact temperature cut, since most of the 
hot gas in groups and cluster has a temperature that is a factor of 


two greater than the virial temperature. This temperature selection 
is used as a rough proxy for the hot X-ray emitting gas. A proper 
selection of X-ray emitting gas is beyond the scope of this paper, 
as this would require an accurate computation of the X-ray lumi¬ 
nosity of the gas particles. A luminosity weighted scheme for the 
shape of the hot gas would result in the inner regions domin ating 
the shape resulting in more spherical shapes JCrain et al.l2013i ). On 
the other hand, a mass weight scheme, as adopted in this work, 
would be closer to the shape that a Sunyaev-Zeldovich (SZ) exper¬ 
iment would measure since SZ flux is proportional to the gas mass 
and the temperature, making it potentially testable with a combined 
SZ-lensing analyses. 

Fig.[7]pr e sents the variation of the shape parameters, spheric¬ 
ity on the left and triaxiality on the right, of the temperature- 
selected hot gas particle. The convergence of the sphericity param¬ 
eter between the different simulations is poorer in this case than 
for other components shown earlier. By imposing a strict limit on 
the number of particles needed for measuring the shape we limit 
our results to only few points for the EAGLE simulations. For in¬ 
stance, is no longer possible to connect the results from L025 and 
L100. Nonetheless by relaxing the constraint on the number of par¬ 
ticles (grey points), it is possible to identify a trend in the shapes 
that suggests an increasing triaxiality and decreasing sphericity of 
the hot gas component with host halo mass. 

We have also studied the radial dependence of the shape pa¬ 
rameters for the hot gas component of haloes (not shown). Given 
the limit on the minimum number of particles, only three mass bins 
could be investigated (M|qo > 10 12 /i _1 Mg) and only down to 
radius of r/Goo — 0.3, for which no significant radial trend was 
found. 


5 MISALIGNMENT OF GALAXIES WITH THEIR OWN 
HOST HALOES 

In this section, we show the relation between the orientation of 
haloes and that of their stellar and hot gas component. Specifically, 
we will show how the orientation of the major axis of the stellar 
distribution (S j5.lt and of the hot gas distribution ( f|5.2b compare to 
that of the host halo. Similarly to the case of the shape parameters, 
we will investigate the mass, radial, and redshift dependence of this 
relation. Note that we focus our study mainly on central haloes. We 
remind the reader that a formal definition of the axes of particle 
distribution and their relative misalignment angles is provided in 

es 

5.1 Misalignment of stars with their host haloes 

The left panel of Fig. [8] shows the median misalignment of stars in 
spheres of increasingly larger radii with the direction of the total 
matter distribution within the virial radius for different bins in halo 
mass and for radii expressed in units of Goo - Perhaps not surpris¬ 
ingly, the alignment of stars within the total halo increases from 
the inner to the outer part of the halo. The gradient is relatively 
steep, with the misalignment angle between the stars and their host 
haloes decreasing from about 30 degrees (at r ~ 0.03r2oo) to a 
few degrees (at r ~ r^oo) in the case of the most massive haloes. In 
less massive haloes, the misalignment is larger at all scales. Similar 
trends hold at 2 = 1 (dotted lines). The right panel of Fig. [8] shows 
the misalignment of stars with the direction of the total (mostly 
dark) matter, where both are now enclosed in spheres of increas¬ 
ingly larger radii. At each radius, the misalignment is small. Stars 
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Figure 6. Stellar shape parameters (sphericity on the left, triaxiality on the right) as a function of distance from the centre of the halo for haloes in different 
halo mass bins (see legend). Distances have been rescaled to the mean halo radius in each mass bin, to ease the comparison of the results for different 
masses. Different colours indicate different simulations, whereas different line styles refer to different redshifts. The distribution becomes less spherical and 
more prolate with increasing distance from the halo centre. Vertical arrows indicate the median values of the half-mass radii in stars, which can be 

considered a proxy for the typical extent of a galaxy. The blue arrow lies beneath the green one. 




Figure 7. Shape parameters (sphericity on the left, triaxiality on the right) of the gas distribution for the hot (T > 10®K) component. Different colours 
indicate different simulations, whereas different line styles refer to different redshifts. The error bars represent one sigma bootstrap error on the median. Grey 
lines show the results for mass bins containing haloes with less than 300 hot gas particles. 


are aligned with the total mass to within a few degrees in the most 
massive haloes, whereas the alignment deteriorates to about 20-30 
(10-20) degrees for the least massive haloes at z = 0 (z = 1). 

The misalignment of stars with their host halo can vary sub¬ 
stantially depending on the radius and the mass of a halo. The ar¬ 
rows in the plot represent the values, in units of rloo, of the half 
mass radius in stars, which is a good indicator of the physical ex¬ 
tent of a galaxy. At this radius the orientation of the galaxies is 
clearly a biased proxy of the orientation of the halo. Galaxies are, 
however, much better aligned with the local distribution of mat¬ 
ter. This indicates that the the stellar orientations follow that of the 


dark matter, which is the dominant component in mass, and the 
dark matter itself changes orientation from the inner to the outer 
halo. This causes the stars to be well aligned with the local mass 
distribution but misaligned with the orientation of the entire halo. 

5.1.1 Probability distribution function of misalignment angles 

In the previous section we presented the median value of misalign¬ 
ment between the halo and the stellar component. The upper panel 
of Fig.|9]shows the probability distribution function of the cosine of 
the misalignment angle between the stars and the entire host halo 
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Figure 8. Spatial variation of the median cosine of the misalignment angle between the major axis of stars and the underlying (mostly dark) matter distribution. 
Different colours indicate different halo mass bins, whereas different line styles indicate different redshifts. Radial coordinates are normalized by the mean 
halo radius, t^qq , of each mass bin to ease the comparison. Only central haloes are used. Vertical arrows represent the median value of rjjjj in units of 
in different mass bins. Left Panel. Median value of the cosine of the angle between the major axes of the stellar component and that of the entire halo. Here 
the direction of the halo is determined using all particles belonging to the halo. Right Panel. Median value of the cosine of the angle between the major axes of 
the stellar component and that of the halo. The misalignment between the stars and halo is caused, to first order, by the misalignment of the inner dark matter 
halo with the total matter distribution in the halo. 
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Figure 9. Upper panel: probability distribution function of the cosine of the misalignment angle between the major axis of the distribution of stars inside 
r h£df ’ a °d the major axis of the entire halo for four halo mass bins. The black histograms indicate the probability distributions for the total sample of haloes 
that satisfies the resolution criteria, whereas coloured histograms refer only to the subsample of haloes whose mass is indicated in the legend. Vertical lines 
indicate the median values of the misalignment angle (same colour convention as for the histograms). Red dashed curves represent the analytic fit discussed in 
Appendix [B] Lower panel: cumulative version of the probability function for early- and late-type galaxies (dotted and dashed curves respectively). 
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Figure 10. Probability distribution function of the 2D misalignment angle between the major axes of the projected distribution of stars (inside rj'j and the 
major axes of the projected total matter distribution for four halo mass bins. The black histograms indicate the probability distributions for the total sample 
of haloes that satisfies the resolution criteria, whereas coloured histograms refer only to the subsample of haloes whose mass is indicated in the legend. 
Vertical lines indicate the median values of the misalignment angle (same colour convention as for the histograms). Red dashed curves represent the analytic 
fit discussed in Appendix[B] whereas black and grey curves are obtained with analytic functional forms that have been employed in the literature (see text). 


for central galaxies. Here the stars are taken to be inside . Each 
panel shows a different mass bin and therefore a different simula¬ 
tion. The colour histograms show the misalignment distribution for 
haloes in that specific mass bin, whereas the black histograms show 
the probability distribution functions for all haloes that are above 
the halo mass resolution limit (300 stellar particles inside rf^f) in 
the corresponding simulation . The vertical lines show the median 
values for the distributions and the dashed red curves are analytic 
fits (see Appendix iBt. The lower panel of Fig. [9] shows the cumu¬ 
lative probability of the cosine of the misalignment angle for early- 
(dotted curves) and late-type galaxiesQ (dashed curves) as well as 
for the whole sample of haloes (continuous curves). 

The distribution of the cosine of the misalignment angle has a 
long tail towards low values (i.e. strong misalignment) with a floor 
value that decreases with increasing halo mass. The misalignment 
angle distribution of resolved haloes is quite similar in shape for the 
different simulations. Using the fitting functions provided in Ap¬ 
pendix [B] and the median values of the misalignment angle shown 
in the previous plots, it is possible to populate dark matter haloes 
with galaxies oriented such that these misalignment distribution are 
reproduced. 

iBett et all d20ld ) quantified the misalignment angle between 
the stellar and total matter distribution in a sample of (about 90) 
disc galaxies selected from a hydrodynamic simulation in a cu¬ 
bic volume of 35 h~ x Mpc by side. They found that half of these 
galaxies have a misalign ment angle larger t h an 45 degrees . Using 
the GIMIC simulations JCrain erai]|2009bl ). iDeason et al] J201 ll) 
reported that 30% of disc galaxies with average halo mass of 
logio(M su b/[M 0 ]) = 12.1 have a misalignment angle of more 
than 45 degrees. Both these studies are in broad agreement with 
our findings for similar halo masses. Specifically, in the EAGLE 
simulations, we find that half of the disc galaxies have misalign¬ 
ment angles larger than 50 (40) degrees in L025 (L100) and 30% 
of the g alaxies in L 100 (for which the typical halo mass is close to 
that in lPeason et al] d201 lb have misalignment angles larger than 
60 degrees. 

Fig.llOlshows the probability distribution function of the mis¬ 
alignment between the major axes of the projected halo and the 


4 See definition of disc galaxies in 15.1.21 


projected stellar mass component. For comparison, we report with 
black (me an) and grey (one sigma de viation) dashed lines the re¬ 
sults from lOkumura. Jing & Lif | |2009|) who found that, by assum¬ 
ing a Gaussian misalignment distribution between LRGs and dark 
matter haloes, they were able to account for the discrepancy be¬ 
tween the measured orientation correlation of LRGs and the one 
predicted by N-body simulations. Furthermore, we overplot ana¬ 
lytic fits to our discrete distributions using a double Gaussian (red 
dashed curves, see Appendix Notably, none of our probabil¬ 
ity distributions resembles a Gaussian function. It is obvious that a 
single Gaussian function cannot be used as a fair description of the 
probability functions measured from our simulations. 

5.1.2 Misalignment for early- and late-type galaxies 

In this section, we study the alignment between stars and their 
host haloes in early- and late-type galaxies. Given the galaxy stel¬ 
lar velocity dispersion, <r s tar, and the halo maximum circular ve¬ 
locity, one can define the ratio p = cr a tar/Vj.™c X to quan¬ 

tify whether a galaxy is supported either by ordered (rotational) 
motion or by the velocity dispersion. We adopt the convention 
that p 0.5 indicates a rotationally-supported galaxy (late type), 
whereas p > 0.5 indicates a dispersion-supported galaxy (early 
type). 

The left panel of Fig. ELD shows the median misalignment 
of the direction of the entire host halo with that of stars in 
spheres of increasingly larger radii for early- (dotted lines) and 
late-type (dashed lines) galaxies. As for the entire galaxy popu¬ 
lation, the misalignment of stars with their host halo decreases 
from the inner to the outer part of the halo. The misalign¬ 
ment decreases with mass and is lower for late- than for early- 
type galaxies. The misalignment of early-type galaxies in low- 
mass haloef] (11 < logio(M 2 oo/[/i 1 M 0 ]) < 12 and 12 < 
logio(AL 200 /[ /t. 1 M@] < 13)) is especially large at all radii and 
its radial dependence is significantly steeper than in all other cases. 

The right panel of Fig. m shows the ratio, 

5 In Fig. ED we use the EAGLE L100 simulation also for the least massive 
bin (cyan lines) to improve the otherwise poor statistics of the EAGLE L025 
simulation. 
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Figure 11. Left Panel. Spatial variation of the cosine of the median misalignment angle between stars and the underlying (mostly dark) matter distribution. 
Different colours refer to different halo mass bins, whereas different line styles refer to different galaxy types. Radii are rescaled to the mean halo radius, r^QQ, 
to ease the comparison of the results corresponding to different halo mass bins. Only central haloes are used in order to remove effects that can alter mostly the 
alignment of satellites. Here the direction of the halo is determined using all particles belonging to the structure. A kinematic classification has been employed 
(see text) to divide galaxies into early- and late-type. Right Panel. Spatial variation of the ratio between the alignment of late and early type galaxies. As in the 
left panel, the alignment is expressed in terms of the cosine of the angle between galaxy and halo major axes. In both panels the least massive bin is taken from 
the L100 simulation to improve on the otherwise poor statistics of L025. The vertical arrows represent the median values of in units of r'ff in different 
mass bins. 


cos 0star/cos of the cosine of the misalignment angle 

between the stars of early- and late-type galaxies and the entire 
halo. At all radii of interest here, early-type galaxies are more 
misaligned than late-type galaxies. The misalignment angle of 
late-type galaxies is smaller by about 10-20% at r ~ 0.03r2oo 
approximately the expected physical extent of the galaxy. 

A more detailed investigation of the galaxy-halo misalignment 
as a function of galaxy type is beyond the scope of this paper. We 
do acknowledge that this is certainly an interesting direction to be 
further explored, especially in view of the fact that many (current 
and forthcoming) lensing studies for which the misalignment angle 
hampers the interpretation of the signal use early-type galaxies such 
as LRGs. This exploratory work suggests that late-type galaxies are 
instead less misaligned with their host halo and therefore, in this 
respect, to be preferred to early-type galaxies. 


5. 1.3 Misalignment for central and satellites galaxies 

In this section, we characterize the alignment of stars with their 
host halo for centrals and satellites separately. Note that we 
have only considered centrals in the preceding sections. The left 
panel of Fig. m shows the median misalignment of the direc¬ 
tion of the entire halo with that of the stars in central and satel¬ 
lite galaxies, whereas the right panel of Fig. ED shows the ra¬ 
tio of the cosine of the misalignment angle between the en¬ 
tire halo and the stars for central and satellite galaxies. As 
for Fig. El we employ the EAGLE LI00 simulation for the 
mass bins (11 < logio(M su b/[/i _1 Mq]) < 12 and 12 < 
logio(M su b/[/i^ 1 Mq]) < 13) to improve the otherwise poor 
statistics of the EAGLE L025 simulations. Furthermore, we adopt 
here the dark matter half-mass radius, rflaTf, as a definition of 
the extent of a halo, as this is properly defined for both cen¬ 


trals and satellites whereas an overdensity with respect to a back¬ 
ground/critical value is an ill-defined concept for subhaloes that 
host satellite galaxies. At all radii, the misalignment angle between 
the entire halo and the stars in central and satellite galaxies is the 
same to within 10%. The radial trend is in qualitative agreement 
with those of the whole sample shown in Fig. [8] (i.e. the misalign¬ 
ment decreases from the inner to the outer halo). 

The consistently lower misalignment in the outer parts of 
satellites could be due to the tidal stripping removing the outer 
(and more misaligned) part of the halo. Instead, in the inner part 
the resulting reduction of TbSf (for which the radii are normalized) 
would produce a shift of the whole relation to the right, effectively 
increasing the misalignment. The competition between these two 
processes could explain the transition between a more misaligned 
inner part to a less misaligned outer part of satellites with respect 
to central subhaloes. 

This result indicates that satellite-specific physical processes 
(e.g. dynamical friction, tidal stripping) generally do not have a 
strong impact on the misalignment between the stellar and (mostly 
dark) matter component. 


5.1.4 The effect of projection on the misalignment angle 

Observationally one only has access to quantities projected onto 
the plane of the sky. Therefore, it is of interest to compute the mis¬ 
alignment in a (random) two-dimensional (2D) plane onto which all 
particles of the simulations have been projected. Correspondingly, 
one has 2D inertia tensors that describe the matter distribution of 
each component. In this 2D application, the misalignment angle be¬ 
tween the stars and the halo is measured as the angle between the 
main eigenvectors of the inertia tensor of stars and (mostly dark) 
matter. 
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Figure 12. Same as Fig. ll 1 I but for central and satellite galaxies and subhaloes (see text). To ease the comparison for results of different halo types, radii have 
been rescaled to the half mass radius for the dark matter mass, rjJjJJf • I n both panels the least massive bin is taken from the LI 00 simulation in order to improve 
the statistic. The vertical arrows represent the median value of rjjjj in units of the different mass bins. 
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Figure 13. Comparison between the spatial variation of the 3D (continuous lines) and 2D (dashed lines) of the cosine of the median misalignment angles of 
stars with the underlying (mostly dark) matter distribution. Different colours refer to different halo mass bins. Only central haloes are used to exclude effects 
that can alter mostly the alignment of satellites. As discussed in the main body of the paper, the projection, by reducing the degrees of freedom of the system, 
increases the alignment. 


The left panel of Fig. |T3] shows the radial- and mass- 
dependence of the median (cosine of the) misalignment angle for 
the 3D (solid) and the 2D (dashed) case. Clearly, the net effect of 
projecting the 3D distribution onto a 2D plane is an increase in the 
alignment at all radii and all halo masses. The right panel of Fig. ED 
shows the ratio between the cosine of the misalignment angle in 2D 
and 3D. The ratio decreases with both mass and radius but is always 
greater than unity. It reaches values of about 1.25-1.35 for the low- 
mass bins at the radii that are representative of the physical extent 
of a galaxy. A similar result was reported in Tenl4. 


5.2 Misalignment of hot gas with its host halo 

Fig.ll4lshows the radial and mass dependence of the alignment of 
the hot component of the gas (T > 10 6 K) with its host halo. The 
results are only shown for three mass bins, because the mass bin 
11 < logio(M 2 oo/[ h,- 1 Mpc]) < 12 does not contain enough 
hot gas particles to retrieve reliable estimates for the orientation. 
For the highest halo mass bins (right panel) the spatial variation of 
the misalignment angle between the hot component and the entire 
halo (left panel) is similar to that of the stars in the same halo mass 
bins. On the other hand, the misalignment between the hot gas and 
the local matter distribution differs from the case of stars: the hot 
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Figure 14. Same as Fig.[8]but for the hot (T > 10®K) gas component of haloes. The alignment of the hot component increases with radius and mass. Except 
for the highest-mass haloes, the gas does not follow the dark matter distribution as well as was the case for the stars (c.f. Fig. [8}. 


gas component is significantly misaligned with respect to the local 
matter distribution. Specifically, for haloes in the mass range 12 < 
logio(M 2 oo/[^ 1 M@]) < 13 the misalignment angle of the hot 
gas is as large as 50 degrees at r ~ 0.3 r^oo and it is ~ 30 (~ 10) 
degrees for the halo mass range 13 < logio(M 2 oo/[/i 1 Mg]) < 
14 (14 < logio(Af 2 oo/[/t _1 M 0 ]) < 15). Results for redshift 
2 = 1 (dotted lines) have similar radial and mass dependence as 
for redshift 2 = 0. 

Because it is observable out to larger radii than the stellar dis¬ 
tribution of the central galaxy, hot gas represents a valuable tracer 
of the gravitational potential of massive clusters. Unfortunately, the 
fact that the hot gas tends to be largely misaligned with the local 
matter distribution makes it a poor tracer of the shape of the halo, 
unless logio(M 20 o/[/i _1 Mg]) > 14. 


6 SUMMARY AND CONCLUSIONS 


This paper reports the results of a systematic study of 
halo an d galaxy shapes and their relative alignment in the 


EAGLE (lSchaye etal 


2015^ Crain et aill2QI5l i and cosmo-OWLS 
dLe Brun et alJ l20l4 [McCarthy et al.1 |20I4|) hydro-cosmological 
simulations. Several aspects of these simulations make them an 
ideal tool for this investigation. First, the combination of these 
simulations allows us to apply our study to four orders of mag¬ 
nitude in halo masses with sufficient resolution and statistics. Sec¬ 
ond, the EAGLE simulations have been calibrated to be in agree¬ 
ment with the observed present-day galaxy stellar mas s function 
and the observed size-mass relation ISchave et al .1120 1 51) . Third, it 
has been shown that cosmo-OWLS simulations reproduce key (X- 
ray and optical) observed properties of galaxy groups as well as 
the observed galaxy mass function for haloes more massive than 
log(M/[/U 1 M 0 ]) = 13. 

We have studied the shapes of the distributions of dark 
matter, stars and hot gas in haloes with masses 11 < 

logio(M 2 oo/[ fiT 1 M 0 ]) < 15 and their evolution in the redshift 
range 0^2^ 1. We find that the matter distribution in haloes is 
more aspherical (and triaxial) at higher halo mass and higher red¬ 


shift (see Fig.0. The same qualitative trends hold for the star and 
the hot gas distribution in haloes (see Figs.[4]and[7]l. We report (in 
Fig. [6]) the spatial variation of the median of the shape parameters 
of the stellar distribution from ~ 0.02r 2 oo (i.e. a few to tens of 
kpc) to r 2 oo (i.e. up to a few Mpc). We note that at fixed radius 
and halo mass, stellar distributions are generally less spherical than 
dark matter haloes. We have measured the r.m.s. of the projected 
stellar ellipticity as a function of halo mass. We find a modest mass 
dependence, with r.m.s. stellar ellipticity increasing by 50 % as halo 
mass increases by four orders of magnitude. We note that the values 
of the r.m.s. stellar ellipticity vary from ~ 0.2 to ~ 0.35 when one 
considers only stars within the star half-mass radius. However, the 
same quantity varies from ~ 0.35 to ~ 0.55 when all stars within 
the halo are considered (see Fig.[5]l. 

iTenneti et ~al] 120141) recently used the Massive Black II simu¬ 
lation to study the mass dependence and evolution of the stellar and 
dark matter components of haloes and subhaloes. Their findings 
are, for the most part, in qualitative agreement with ours. However 
we find a few differences as reported in the corresponding sections 
(see e.g. § [4] and the discussion of Fig. 0. Specifically, we high¬ 
lighted sources of potential biases in their analysis. As detailed in 
§ [3] those biases mostly stem from the use of a hydro-simulation 
that does not reproduce the observed stellar-halo mass relation and 
by imposing an artificial cut-off in the minimum stellar mass for 
which the shape is calculated. 

We have measured the misalignment of the baryonic compo¬ 
nents (stars and hot gas) of galaxies with their own host haloes. We 
find that stars align well with the underlying (mostly dark) matter 
distribution, especially when all stars inside the halo are considered 
(see Fig. 0. However, the stellar distributions in the inner parts of 
the host haloes do exhibit a median misalignment of about 45-50 
degrees. The misalignment is smaller in more massive haloes (13 ^ 
logio(M 2 oo/[ tT 1 M 0 ]) < 15), late-type galaxies (see Fig.lUll. 
and central galaxies (see Fig. EDl. The hot gas distribution can only 
be traced with a sufficient number of particles only in the outer part 
(^ O.3r 200 ) of massive (12 ^ logio(M 20 o/[/i _1 M 0 ]) ^ 15) 
haloes. In this range we find that the alignment of the hot gas with 
the entire halo is similar to that of the stellar distribution. However, 
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the hot gas does not align well with the local matter distribution, ex¬ 
hibiting misalignment angles larger than 20 (typically 30 to 50) de¬ 
grees in haloes with masses 13 ^ logio(M 2 oo/[/i _1 Mq]) SC 15 
(see Fig.mil 

We have quantified the effect of projection on the median mis¬ 
alignment angles between the stellar distribution and the halo (see 
Fig.ED. Projection reduces the degrees of freedom of the system, 
increasing the alignment. Finally, we provided the probability dis¬ 
tribution of the misalignment angle between the major axis of the 
stellar distribution inside the stellar half-mass radius and the major 
axis of the entire halo for the three- and two-dimensional case (see 
Figs.[9]and[T0] respectively). 

We have encapsulated our results in fitting functions (see 
Appendix B) and tables that allow interested practitioners to 
straightforwardly include our results into halo catalogues ex¬ 
tracted from N-body simulations. The complete list of fit¬ 
ting parameters as well as tabulated values are available at 
http: //www. strw. leidenuni v. nl/M V 15 a/ 

A natural extension of this work is the study of the correlation 
functions of galaxy shapes. We will present such an investigation 
in a future publication. 


ACKNOWLEDGEMENTS 

We thank the anonymous referee for insightful comments that 
helped improve the manuscript. MV and MC thank Flenk Hoek- 
stra and Rachel Mandelbaum for useful and stimulating discus¬ 
sions. MC acknowledges support from NWO VIDI grant num¬ 
ber 639.042.814 and ERC FP7 grant 278594. RAC is a Royal 
Society University Research Fellow. This work used the DiRAC 
Data Centric system at Durham University, operated by the Insti¬ 
tute for Computational Cosmology on behalf of the STFC DiRAC 
HPC Facility (www.dirac.ac.uk). This equipment was funded by 
BIS National E-infrastructure capital grant ST/K00042X/1, STFC 
capital grant ST/H008519/1, and STFC DiRAC Operations grant 
ST/K003267/1 and Durham University. DiRAC is part of the Na¬ 
tional E-Infrastructure. We also gratefully acknowledge PRACE for 
awarding us access to the resource Curie based in France at Tres 
Grand Centre de Calcul. This work was sponsored by the Dutch 
National Computing Facilities Foundation (NCF) for the use of 
supercomputer facilities, with financial support from the Nether¬ 
lands Organization for Scientific Research (NWO). The research 
was supported in part by the European Research Council under 
the European Union’s Seventh Framework Programme (FP7/2007- 
2013) / ERC Grant agreements 278594-GasAroundGalaxies, GA 
267291 Cosmiway, and 321334 dustygal, the UK Science and 
Technology Facilities Council (grant numbers ST/F001166/1 and 
ST/1000976/1), Rolling and Consolodating Grants to the ICC, 
Marie Curie Reintegration Grant PERG06-GA-2009-256573. TT 
acknowledge the Interuniversity Attraction Poles Programme initi¬ 
ated by the Belgian Science Policy Office ([AP P7/08 CHARM]) 


REFERENCES 

Bailin J., Steinmetz M., 2005, ApJ, 627, 647 

Behroozi P. S„ Wechsler R. H„ Conroy C., 2013, ApJ, 770, 57 

Bett P„ 2012, MNRAS, 420, 3303 

Bett P., Eke V., Frenk C. S., Jenkins A., Okamoto T., 2010, MN¬ 
RAS, 404, 1137 

Booth C. M„ Schaye J„ 2009, MNRAS, 398, 53 


Boylan-Kolchin M„ Ma C.-P, Quataert E„ 2006, MNRAS, 369, 
1081 

Bryan S. E., Kay S. T., Duffy A. R., Schaye J., Dalla Vecchia C., 
Booth C. M„ 2013, MNRAS, 429, 3316 
Cen R„ 2014, ApJ, 785, L15 

Chen D. N„ Jing Y. P„ Yoshikaw K„ 2003, ApJ, 597, 35 
Cole S„ Lacey C., 1996, MNRAS, 281, 716 
Crain R. A., McCarthy I. G., Schaye J., Frenk C. S., Theuns T., 
2010, ArXiv e-prints (arXiv: 1011.1906) 

Crain R. A., McCarthy I. G., Schaye J., Theuns T., Frenk C. S., 
2013, MNRAS, 432, 3005 
Crain R. A. et al., 2015, MNRAS, 450, 1937 
—, 2009a, MNRAS, 399, 1773 
—, 2009b, MNRAS, 399, 1773 

Croft R. A. C., Di Matteo T., Springel V., Hernquist L., 2009, 
MNRAS, 400, 43 

Dalla Vecchia C., Schaye J., 2008, MNRAS, 387, 1431 
—, 2012, MNRAS, 426, 140 

Davis M„ Efstathiou G„ Frenk C. S„ White S. D. M„ 1985, ApJ, 
292, 371 

Deason A. J. et al„ 2011, MNRAS, 415, 2607 
Dolag K., Borgani S., Murante G., Springel V., 2009, MNRAS, 
399, 497 

Dubinski J., 1998, ApJ, 502, 141 

Hahn 0., Teyssier R„ Carollo C. M„ 2010, MNRAS, 405, 274 
Heymans C., Brown M., Heavens A., Meisenheimer K.. Taylor 
A., Wolf C., 2004, MNRAS, 347, 895 
Heymans C. et al., 2013, MNRAS, 432, 2433 
Hoekstra H., Franx M., Kuijken K., 2000, ApJ, 532, 88 
Jing Y. P„ Suto Y„ 2002, ApJ, 574, 538 

Joachimi B., Semboloni E., Bett P. E., Hartlap J., Hilbert S., Hoek¬ 
stra H„ Schneider P, Schrabback T„ 2013, MNRAS, 431, 477 
Kaiser N„ 1984, ApJ, 284, L9 

Kang X., van den Bosch F. C„ Yang X., Mao S„ Mo H. J„ Li C., 
Jing Y. P„ 2007, MNRAS, 378, 1531 
Komatsu E. et al., 2011, ApJS, 192, 18 

Le Brun A. M. C., McCarthy I. G., Schaye J., Ponman T. J., 2014, 
MNRAS, 441, 1270 

Mandelbaum R. et al., 2014, ApJS, 212, 5 
McCarthy I. G„ Le Brun A. M. C., Schaye J., Holder G. P, 2014, 
MNRAS, 440, 3645 

McCarthy I. G. et al., 2010, MNRAS, 406, 822 
Miller L. et al., 2013, MNRAS, 429, 2858 
Moster B. P„ Naab T„ White S. D. M„ 2013, MNRAS, 428, 3121 
Naab T.. Khochfar S„ Burkert A., 2006, ApJ, 636, L81 
Navarro J. F„ Frenk C. S., White S. D. M„ 1997, ApJ, 490, 493 
Okumura T„ Jing Y. P„ Li C., 2009, ApJ, 694, 214 
Planck Collaboration et al., 2014, A&A, 571, A16 
Reyes R., Mandelbaum R., Gunn J. E., Nakajima R., Seljak U., 
Hirata C. M„ 2012, MNRAS, 425, 2610 
Romano-Dlaz E., Shlosman I., Heller C., Hoffman Y., 2009, ApJ, 
702, 1250 

Rosas-Guevara Y. M. et al., 2013, ArXiv e-prints 
(arXiv:1312.0598) 

Sales L. V., Navarro J. F., Theuns T., Schaye J., White S. D. M., 
Frenk C. S„ Crain R. A., Dalla Vecchia C., 2012, MNRAS, 423, 
1544 

Scannapieco C., White S. D. M., Springel V., Tissera P. B., 2009, 
MNRAS, 396, 696 

Schaye J. et al., 2015, MNRAS, 446, 521 

Schaye J., Dalla Vecchia C., 2008, MNRAS, 383, 1210 

Schaye J. et al., 2010, MNRAS, 402, 1536 


© 2014 RAS, MNRAS OOO.mflsl 


Galaxy-halo misalignment 17 


Sharma S., Steinmetz M., 2005, ApJ, 628, 21 
Springel V., 2005, MNRAS, 364, 1105 

Springel V., White S. D. M., Tormen G., Kauffmann G., 2001, 
MNRAS, 328, 726 

Tenneti A., Mandelbaum R., Di Matteo T., Feng Y., Khandai N., 
2014, MNRAS, 441,470 

Tenneti A., Singh S., Mandelbaum R., Matteo T. D., Feng Y., 
Khandai N„ 2015, MNRAS, 448, 3522 
van den Bosch F. C., Abel T., Croft R. A. C., Hernquist L., White 
S. D. M„ 2002, ApJ, 576, 21 

Viola M„ Kitching T. D„ Joachimi B„ 2014, MNRAS, 439, 1909 
Wang Y., Yang X., Mo H. J., Li C., van den Bosch F. C., Fan Z., 
Chen X., 2008, MNRAS, 385, 1511 
Wiersma R. P. C., Schaye J., Smith B. D„ 2009, MNRAS, 393, 99 
Wiersma R. P. C., Schaye J., Theuns T., Dalla Vecchia C., Toma- 
tore L„ 2009, MNRAS, 399, 574 
Zemp M., Gnedin O. Y., Gnedin N. Y., Kravtsov A. V., 2011, 
ApJS, 197, 30 


APPENDIX A: CAVEATS IN SHAPE PARAMETER 
ESTIMATION 

A1 The choice of inertia tensor 

There exists a plethora of methods designed to characterize the 
shape of a given three-dimensional particle distribution (say dark 
matter, star, gas) in the context of cosmological structure forma¬ 
tion simulations (see Zemp et al. 2011 and references therein). All 
those methods are based on the idea that structures can be well 
described by an ellipsoidal shape. However, the actual algorithms 
used to retrieve this shape can differ substantially and unfortunately 
the corresponding results do not often agree (see Zemp et al. 2011 
for an analysis of this problem under controlled conditions with 
known shapes). Most notably, results on the shape of a particle dis¬ 
tribution may vary if one adopts the inertia tensor rather than the 
reduced in ertia tensor, or some iterative form of the two (see dis¬ 
cussions in lZemp et alJl20lTI : lTenneti et alJl2015il . The differences 
between the inertia tensor and the reduced inertia tensor are driven 
by the fact that in the reduced inertia tensor calculation particles are 
not weighted by their distance from the centre. The net effect is that 
if the reduced inertia tensor is used, the shape is less dominated by 
the particles in the outer part of haloes, meaning that the retrieved 
shape tends to be more spherical as particles in the inner parts of 
haloes are typically more spherically distributed. We repeated our 
analysis using the reduced inertia tensor and found that there is lit¬ 
tle information content in exploring the radial variation using this 
method since its properties have almost no variation with radius. 
We also note that the misalignment angle, which is the quantity of 
primary interest here, is less affected than shape by the choice of the 
algorithm that defines the shape parameters. This is especially true 
when the alignment is calculated between particles distributions at 
the same distance from the centre. 

Another possible variation in the shape calculation is to use an 
iterative method for the inertia tensor calculation but this method 
was proven to gi ve very similar res ults when the inertia tensor is 
used, as shown in lTenneti et al.112015l) . 

Throughout the paper we adopted the definition presented in 
eq. 0 and the corresponding shape parameters defined in eq. 0 . 
While our adopted method may be considered somewhat arbitrar y 
(e.g. see the discussions in I.Tine & Sutoll2002l : IZemn et alJl201 lh . 


we used this a pproach as it is adequate for the comparison we pre¬ 
sented (see e.e lBettl2012l) and because it allowed us to compare our 
results with most of the other results in the literature. 


A2 The effect of sampling 

An important technical aspect regarding shape measurements is to 
find the minimum number of particles required to obtain a reliable 
estimate. To this aim, we simulate a three-dimensional halo with a 
given axis ratio and a Navarro, Frenk, & White (hereafter NFW) 
density profile l lNavarro. Frenk & Whitej|T997t) . More specifically, 
we choose values for the three-dimensional halo axis a, b and c, and 
use an analyitical NFW profile with c = 5 and r v i r = a (the largest 
axis). We then generate A part spherical coordinates (r, <fi, 6) using 
the NFW profile as a selection function, redrawing any coordinates 
that fall outside the ellipsoid defined by a, b, and c. Specifically, 
we use 1 ^ Apart ^ 3000. For each value of Apart, we repeat 
the sampling 10 5 times so as to obtain a median and a standard 
deviation. 

It is worth noting that the number of particles needed for an 
unbiased shape determination depend on the intrinsic shape of the 
halo. Many more particles are needed to retrieve a quasi-spherical 
shape than for example a disky structure. For our test, the intrin¬ 
sic shape of the halo was chosen to have sphericity S = 0.6 and 
triaxiality T = 0.7, which is representative of the average shape 
parameters of our halo sample (see e.g. results in i ]4.2.2l and Fig.0. 

In Fig. IA1I we show the relative error on the retrieved shape 
parameters, S (green lines) and T (red lines), as a function of 
A par t. Solid lines refer to the median, whereas dashed lines refer 
to the 16th and 84th percentiles. The retrieved sphericity shows 
a monotonic trend with the number of test particles. The spheric¬ 
ity increases towards the real value as the number of test particles 
is increased. This means that any resolution effect will lead to an 
underestimating of the true sphericity of haloes. For this particular 
halo shape using 300 particles will lead to an average ~ 2% error 
in the determination of the sphericity with an accuracy of ~ 10%. 
The triaxiality is typically underestimated but converges faster to 
the true value, with the systematic error dropping below 3% for 30 
particles. On the other hand the scatter around the median converg¬ 
ing slowly and is still 20% for 300 particles. Triaxiality thus re¬ 
quires more particles than sphericity in order to reduce the random 
error below a specific value. Throughout the paper, we thus employ 
A par ts ^ 300 as the limit for shape parameter determination. This 
assures very good estimate of the median value of the shape pa¬ 
rameters with a systematic error below 3% and a random error of 
10% in the sphericity and 20% in triaxiality. In this work we did 
not show these systematic errors in the shape measurement but only 
the statistical errors evaluated using the bootstrapping technique. 


APPENDIX B: ANALYTIC FITS FOR THE 
MISALIGNMENT ANGLE DISTRIBUTIONS 

In this section we provide fitting function^ for the distribution of 
the cosine of the 3D misalignment angle 9, as well as for the 2D 
misalignment angle, 0. We note that the choice of using the cosine 
as the variable of the fitting function stems from the notion that the 

6 The analytic fits provided in this section reproduce the median of the 
distributions obtained from the simulations with an accuracy better than 

1 %. 
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Figure Al. Convergence test for shape parameter retrieval. Relative error 
on the retrieved shape parameter of a synthetic NFW halo as a function of 
the number of particles used to sample the underlying distribution. The rel¬ 
ative errors on sphericity and triaxiality are indicated by the green and red 
lines, respectively. The test is performed using a typical sphericity value for 
the synthetic halo, S = 0.6 and T = 0.7. For each number of particles, 
distributions are drawn 10 5 times and we report the 50th (continous lines), 
16th and 84th percentiles (dashed lines). Retrieving both shape parameters 
with a systematic error smaller than a few percent requires at least 300 par¬ 
ticles. 

distribution of the cosine of the alignment angle of a random set of 
3D vectors is flat, whereas the distribution of the angle itself is not, 
as it is skewed towards large alignments. 

We employ the following functional form: 

Ai 3D (x) = A + exp [B — A(1 — , (Bl) 

where x = cos (6) and 0 < 8 < ■n/2. This functional form has 
four free parameters: A, B, A, fi. We find this number of param¬ 
eters necessary to adequately reproduce the main features of the 
results obtained from the simulations. In the main body of the pa¬ 
per (see Fig. [9}, we have employed this fitting function to describe 
the misalignment angle between the stellar component and its host 
halo in four halo mass bins and for the typical extent of a galaxy, 
the half mass radius r^Jf. The corresponding fitting parameters are 
given in Table [bTI Parameters that refer to other components, ra¬ 
dius definitions, and halo mass bins, as well as tabulated median 
values, can be found at http://www.strw.leidenuniv.nl/MV15a/ 

We analytically describe the probability function of the cosine 
of the 2D misalignment angle with the following functional form: 

M 2D (x) = Cexp (- 7 ^ 2 ^ + £*exp + E, (B2) 

where C, cri, D, 02 , E are the 5 free parameters required to de¬ 
scribe a double Gaussian plus a ‘floor’. The level of complexity of 
this functional form is motivated by the results obtained from the 
simulations. In the main body of the text (see especially Figl 10b. 
we describe the probability distribution of the 2D misalignment 
angle between stars and their host haloes in four halo mass bins 
and accounting only for stars within the typical extent of a galaxy, 
the half mass radius r^{ { . The corresponding fitting parameters are 
given in Table lB2l Parameters that refer to other components, ra¬ 
dius definitions, and halo mass bins, as well as tabulated median 
values, can be found at http://www.strw.leidenuniv.nl/MV15a/ It 
is instructive to compare these 2D misalignment angle distribu¬ 
tions to the commonly assumed single-Gaussian distribution (see 


Simulation 

mass bin 

A 


B 

A 

P 

EA L025 

[ 11 - 

12] 

1.52 E - 

02 

-3.58 

5.92 

1.01 

EAL100 

[12- 

13] 

6.43 E - 

03 

-0.05 

5.13 

0.15 

CO L200 

[13- 

14] 

4.461? - 

03 

-1.04 

5.64 

0.41 

CO L400 

[14- 

15] 

4.13E - 

03 

-0.53 

7.13 

0.42 


Table Bl. Fit parameters for Eg. IB I l iliat describes the misalignment angle 
distribution between the direction of the stellar component inside and 
that of the entire halo. 


mass 

bin 

cr 1 

<7 2 

C 


D 


E 


[11- 

-12] 

5.00 

28.17 

4.691? - 

02 

4.691? - 

02 

4.69 E - 

- 02 

[12- 

-13] 

5.00 

31.65 

1.31 E - 

02 

1.311? — 

02 

1.311? - 

- 02 

[13- 

-14] 

14.70 

32.52 

3.61 E - 

02 

3.61 E- 

02 

3.61 E - 

- 02 

[14- 

-15] 

9.42 

25.71 

4.281? - 

02 

4.28 E- 

02 

4.28 E - 

- 02 


Table B2. Fit parameters for the double Gaussian fitting function Eq. IbII 
that describes the misalignment angle distribution between the direction of 
the projected stellar component inside /'[''[[[ and that of the entire (pro¬ 
jected) halo. 

e.g. lOkumura. Jing & Li|[2009l) . None of the distributions found in 
this study resembles a single-Gaussian and we therefore caution 
interested practitioners against adopting this assumption. 


APPENDIX C: RESOLUTION TEST 

In this section we make use of our different simulations to test the 
influence of resolution on our results. For this test we make use of 
the fact that the L025 simulation is the high-resolution version of 
L100 simulated using a smaller box size. The same is true for L200 
and L400. We do not compare results from simulations that were 
not run with the same code. 

In Fig. 1C II we show in the upper panels the variation of the 
sphericity of the stellar component. In the left panel we show the 
two mass bins for which is it possible to obtain results for both 
the L025 and L100 EAGLE simulations. On the right we do the 
same for the L200 and L400 cosmo-OWLS simulations. Different 
colours refer to different mass bins where as different line styles 
refer to different simulations. In the lower panels we show in the 
same manner the misalignment between the stellar component and 
the whole halo. 

The convergence is generally good, especially at larger radii, 
even though the box size, and hence the halo samples, also change 
between the different simulations. The only case that shows a rel¬ 
atively poor convergence is the misalignment for the least mas¬ 
sive bin of the cosmo-OWLS simulations (blue lines) for which the 
shape of the curves are similar but the values are shifted between 
the two simulations. 
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Figure Cl. Resolution test for the variation of the sphericity (upper panels) and misalignment with the halo (lower panels) of the stellar component as a 
function of radius. We show separately the results for EAGLE (on the left) and cosmo-OWLS (on the right). For each set of simulations we show the results in 
two distinct mass bins. 
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